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Non-linear Effects in Resonant Tunnelling Through a Quantum Dot 
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We study resonant tunnelling transport properties of a quantum dot connected to two leads. The 
flowing electrons are supposed to interact strongly at the quantum dot. The system is represented by 
an Anderson impurity Hamiltonian. The transport properties are characterised by the well-known 
Coulomb blockade properties. However, superimposed to them, this ID system possesses novel 
non-linear current bistability behaviour and other instability phenomena, which leads eventually to 
chaos. 
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INTRODUCTION 



Transport properties have been extensively studied due 
to their interesting device application and very rich phe- 
nomenology Since the first observation of resonant tun- 
nelling in. a 3D double barrier structure by Chang, Esaki 
and TsuEJ, these devices have been found to present non- 
linear properties, which are reflected in the observation of 
multistabilitics in the I-V characteristic curve. Besides a 
peak due to simple resonant transmission, dynamical in- 
trinsic bistability and hysteresis in the negative differen- 
tial resistance region has been measuredau. These prop- 
erties have been explained as a non-linear effect due 
the Coulomb interactions among the flowing charge; 
The phenomenon can be thought to be produced by the 
rapid leakage of the electronic charge accumulated at the 
well between the barriers when the applied potential is 
just taking the device out of resonanceu. Very recently 
we were able to show that the addition of a magnetic field 
parallel to the current induces self-sustained intrinsic cur- 
rent-oscillations in an asymmetric double barrier struc- 
tured. These oscillations were attributed to the enhance- 
ment due to the external magnetic field of the non-linear 
dynamic coupling of the current to the charge trapped 
in the well. These results showed that the system bifur- 
cates as the field is increased, and may transit to chaos 
at large enough fields. Although the presence of the ex- 
ternal magnetic field introduces some ID character to 
these systems, all the properties are obtained from 3D 
systems. In recent years, electron channels have been 
fabricated with a lateral quantum confine energy levels 
sufficiently separated one from the other as to be able 
to study properties of a ID system. When a quantum 
dot is connected through these electron channels to a 
battery, the system exhibits pronounced periodic oscil- 
lations of the conductance when the state of charge of 
the dot is modified by changing the gate potential. This 
is a consequence of the Coulomb blockade of an elec- 
tron when it tries to go into a quantum island, which 



is already occupied by another electronBQ. Although the 
multistability phenomena found in standard 3D resonant 
tunnelling devices have not been detected in ID systems 
where the Coulomb blockade is dominant, it is evident 
that both phenomena derive from the same Coulomb in- 
teraction acting between the flowing charges. In this pa- 
per we developed a self-consistent theory of resonant tun- 
nelling in quasi-lD systems which predicts that, besides 
the normal Coulomb blockade behaviour in the region 
of negative differential resistance, the system develops 
multistable properties. Moreover we show that, while 
the phenomena of hysteresis persists, the system bifur- 
cates and is capable of further bifurcations as the gate 
potential applied to the dot is increased, leading even- 
tually to true chaos. This fact emphasizes the absence 
of stationary solutions for the current in certain regions 
of the parameter space, indicating the existence of time- 
dependent phenomena, probably oscillations. 



THE MODEL 

A theoretical understanding of these properties should 
consider that the system is in a non-equilibrium situation 
under the effect of many-body interactions. A ID An- 
derson Hamiltonian microscopically describes the system 
where the quantum dot plays the role of the impurity. 
The leads are represented by a ID tight-binding Hamil- 
tonian connecting the dot to two reservoirs characterised 
by Fermi levels £i and £ r . The difference £i — £ r cor- 
responds to the potential drop from left to right along 
the sample. The nearest-neighbour Hamiltonian can be 
written as, 
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where V p is the external gate which controls the energy of 
the localised state within the well, U is the local electronic 
repulsion of the charges inside the well localized at site 
0, t is the nearest-neighbour hopping matrix elements, t r 
and ti are the coupling matrix elements between the well 
and the right and left contacts. It is assumed the sys- 
tem has N sites. The electron-electron interaction will 
be treated within the context of the Hubbard approxima- 
tions, which is adequate to treat the Coulomb blockade 
and the non-linear effects derived from the Coulomb in- 
teraction. However, it will not be capable of describing 
the low-lying excitations associated to the Kondo effect 
which in principle the system would have at tempera- 
tures below the Kondo temperature. On a tight-binding 
basis, the stationary state of energy £k can be written as 
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where 4>i a is a Wannier state localized at site i of spin 
a and the coefficients af a obey the non-linear difference 
equations, 
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[3,4] can be obtained through an adequate iteration of it 
from right to left. For a given transmitted amplitude, the 
associated reflected and incident amplitudes may be de- 
termined by matching the iterated function to the proper 
plane wave at the far left, equation [6]. The transmission 
square modulus of T, adequately normalised by the inci- 
dent amplitude I, obtained from the iterative procedure, 
multiplied by the wave vector k, gives us the contribution 
of this wave vector to the current. Due to the presence of 
the non- linear term, the equations [3,4] have to be solved 
self-consistently. With this purpose, we define a second 
pscudo-time-likc iteration in the following way. Initially 
Eqs. [3,4] are solved ignoring the non-linear term and for 
energies up to the Fermi energy. The coefficients thus ob- 
tained correspond to a solution for non- interacting elec- 
trons. They are used to construct the non-linear term for 
the next solution. The procedure is continued, using for 
the non-linear term solutions corresponding to the pre- 
vious iteration. This procedure defines a spatial-pseudo- 
tcmporal map that is linear in space and non-linear in 
time. Maps non-linear in space and time have been stud- 
ied in other contexts. In particular, the coupled map 
model for open flow has been found to exhibit spatial 
; shads' with temporal periodicityQ. A similar behaviour 
has been found studying resonant tunnelling properties 
of a 3D double barrier devices under the effect of an ex- 
ternal magnetic field along the current (4). As we shall 
show, in our case chaos develops in the iteration of time, 
as the self-consistent loop is developed. 
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NUMERICAL RESULTS 



where n~ = 1— < n 0)Cr > and < n 0)CT > is the number 
of electrons with spin a at the well site and is calculated 
from the tight-binding amplitudes at the well, according 
to the equation, 
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where the sum over k covers all occupied electrons states 
of the system with energy below the Fermi level, incident 
from the emitter side. In order to study the solutions of 
equation [2] we assume a plane wave incident from the left 
with an intensity /, with a partial reflection amplitude 
R. The waveform at the far right is a simple plane wave 
with intensity given by the transmission probability T. 
Taking this to be the solution at the far right and left 
edges of the system, we can write, 
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where L is the length of the active part of the system in 
units of the lattice parameter. The solution of equations 



We study a model which consists of leads connected 
to a quantum dot represented by t r — O.li, ti = O.lt, 
U = l.Oi and a Fermi level e/ — 0.05t. These parame- 
ters could adequately represent a semiconductor device 
constructed from a GaAs substrate. The normalisation 
of the wave function is taken so that each site could have 
up to a maximum of two electrons. The map is linear in 
space so that the normalisation involves a simple multi- 
plication of all coefficients by a constant. Since the elec- 
tron density has rather long-range oscillations, we made 
sure the sample was long enough to make finite size ef- 
fects negligible. This is guaranteed by taking a sample 
of N > 100 sites. 

Fig 1 shows the current voltage (I-V) characteristics 
for several values of the gate potential V p applied at the 
quantum dot. As it is shown in the figure for this ID sys- 
tem, we obtain a similar behavior to the very well-known 
bistability in 3D systems. There are two solutions for 
the current as the voltage is increased when the dot is 
charged (solid line), and as the voltage is decreased when 
the dot is without charge (stars). In either case, the self- 
consistent solution converges to a stable fixed point after 
some finite number of iterations. Similarly to 3D tun- 
nelling, the choice of an asymmetric structure (t r > tf) 
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and an increase of the Fermi level augments the charge 
content of the well and consequently the non-linearities 
and in particular the width of the hysteresis loop for this 
ID system. When the gate potential is greater than a 
threshold value, which depends upon the parameters of 
the system, a completely novel feature starts to develop 
where two periods and two fixed points are encountered 
for certain voltages as shown in Fig lb. This bubble-like 
solution is obtained as the voltage is increased outside 
the region of bistability showing that although both phe- 
nomena are derived from the non-linearities introduced 
by the local Coulomb repulsion, they correspond to two 
different effects. Note that these new solutions are not 
reached, as the voltage is decreased (stars) as it corre- 
sponds to a discharged dot where the non-linearities are 
negligible. As the potential gate is raised still further, 
the area of the bistable bubble increases and it enters 
into the bistable region. At larger gate potentials the so- 
lution undergoes further bifurcations and finally a bifur- 
cation cascade leading to a chaotic region as it is shown 
in figlc and Id for a gate potential V p = —0.7t and — 1.4f 
respectively. 

As the external applied potential increases still further, 
the chaotic region eventually unfolds to a solution of finite 
periodicity or to a unique self-consistent one. In fig. 2 
a phase diagram in the parameter space V — V p shows 
the various structures we have obtained. The thick lines 
define the region of bistability, while the stars delimit the 
region where bifurcations of various orders and chaotic 
behaviour take place. 
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FIG. 1. I — V characteristic curve for a system with 
U = It; t r = U = O.lt for four different values of the gate 
potential; a) V p = -O.lt, b) V p = -0.2t, c) V p = -0.7t, d) 
V v = -1.4t. 



FIG. 2. Phase diagram in the parameter space V — V p 
for the system of Figl. The thick lines and stars delimit the 
region of bistability and instability respectively. 



CONCLUSIONS 



We have shown that a quantum dot connected to two 
leads and under the effect on an external applied poten- 
tial have, superimposed to the Coulomb blockade oscil- 
lations of the current, bistable solutions and instabilities 
which eventually become chaotic. Although bistabilitics 
in the current are a well-known phenomenon in 3D reso- 
nant tunneling, to the best of our knowledge, neither its 
existance has been proposed nor has it been measured 
in a ID system. Bistabilities and general instabilities, 
although arising from the same Coulomb interaction, ap- 
pear in different regions of the parameter space reflecting 
the fact that they are different phenomena. A full-time 
dependent study of this problem is under way. 



ACKNOWLEDGMENTS 



This work was partially supported by the Brazilian 
Agencies CNPQ and Fincp, DGICT-UCN, FONDECYT 
Grants 1980225 and 1960417 (Chile). 



3 




V p (units oft) 



